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Abstract 

We study the vicinal surface with adsorption below the roughening temper- 
ature, using the restricted solid-on-solid model coupled with the Ising model. 
By the product-wavefunction renormalization group method, we calculate the 
surface gradient p and the adsorption coverage G as a function of the Andreev 
field 7] which makes surface tilt. Combining Monte Carlo calculations, we show 
that there emerges effective attraction between the steps. This attractive in- 
teraction leads to instability against step bunching. 

Keywords Equilibrium thermodynamics and statistical mechanics, Ising 
models, Monte Carlo simulations, Faceting, Step formation and bunching, 
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I. INTRODUCTION 



Below the roughening temperature Tr, the vicinal surface which is a slightly tilted surface 
relative to one of the facet plane of a crystal, is well described in terms of terraces, steps 
and kinks (TSK picture). Under the picture, the vicinal surface is known to show Gruber- 
Mullins-Pokrovsky-Talapov (GMPT) universal behavior PJ. The surface free energy (per 
projected area), then, is written as |TTj— p|| 

f(p) = f(0)+lP + Bp 3 + O(p i ). (1.1) 

where p is the step density, 7 is the step tension, and B represents the effect of step-step 
interactions. The step density p relates to the surface gradient p (along a crystal axis) as 
p = tan0 = pa,i {<P is the tilt angle of the vicinal surface, and is the height of a single 
step). For convenience we take the units where = 1 in the present article. To be precise, 
7 and B depend on the mean running direction angle (relative to one of the crystal axes on 
the facet plane) of the steps, which we denote by 9; we should then write 7 = 7(6*, T) and 
B = B(9,T) (T: temperature). 

From the thermodynamics of equilibrium crystal shape, the GMPT form of the vicinal 
surface free energy leads to 

f) = ^=7 + 3Bp 2 + 0(p 3 ), (1.2) 

where rj is Andreev field @ which makes the surface tilt. 

It has been known that adsorbed atoms often change properties of the surface |6|,|8Hi7 



To discuss the adsorption effect on the vicinal surface, we take the restricted solid-on-solid 
(RSOS) model |L8j on the square lattice, coupled with the Ising spin system representing 



the adsorbed gas |L9| . The effect of vapor pressure of adsorption in environmental gas phase 



is taken into consideration through the surface chemical potential of adsorbates. 
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II. THE RSOS-ISING COUPLED MODEL 



In the RSOS model, we restrict each nearest-neighbor (nn) height difference Ah to be 
Ah = 0, ±1, which is a reasonable simplification because configurations with large \Ah\ are 
energetically unfavorable. 

We assume that the gas atoms are likely to adsorb at step edge positions and that 
the adsorbed atom modifies the ledge energy locally. Ferromagnetic interactions (attractive 
interaction, in the lattice-gas picture) with coupling constant J between the nearest-neighbor 
spins on the rotated square lattice is considered. We assume simple linear modification of 
the ledge energy as e — > e(l — aa) (er: Ising spin). This modification leads to the interaction 
between the RSOS model and the Ising model. 

The Hamiltonian of the RSOS-Ising coupled model is, therefore, written as 

H= acr b(i,j))\hi - hj\ 

<i,j> 

-J Yl Wv-H £ a b , (2.1) 

<b,b>> <b,b>> 

where hi is the integer surface height at site i, e the "bare" ledge energy, db(i,j) = ±1 the 
Ising spin variable on the bond b(i,j) connecting the nn site pair < i,j >, and H is an 
external field corresponding to the surface chemical potential for adsorbate. We should note 
that the RSOS condition (\Ah\ < 1 for each nn site pair) is implicit in fl2.1|) . (Fig. 1) 

We set Cartesian coordinate system to the square lattice so that the mean running 
direction of steps becomes y direction. The Andreev field rj []7j along the x-direction to tilt 
the surface is introduced by adding a term — r]Ylm,n(h(m+i,n) ~ fym,n)) m the Hamiltonian 
( P-IP ((m, n) is the position vector of the lattice site). 



The Andreev surface free energy f(rj) which directly gives the equilibrium crystal shape 
(ECS) is calculated from the model Hamiltonian with the Andreev-field term by 



/fo) = — AfcTlnZ, (2.2) 



where Z is the partition function of the RSOS-Ising system. 



We calculated the partition function by the transfer-matrix method. We extend the 



mapping between the RSOS model and the vertex model on the dual lattice, [18|J2(J to 
obtain a "decorated" vertex model (Fig.l) which is regarded as a 19 x 16 = 304 vertex 
model. For approximate diagonalization of the transfer matrix, we employ the product- 
wave-function renormalization group algorithm (PWFRG) pi| , |2"2" | which is a variant of the 



White's density-matrix renormalization group (DMRG) method [^3 . 

In the vertex-model representation, the surface gradient p along the re-direction is just 
the thermal average of the vertical edge variable of the vertex model, which can easily be 
calculated from the fixed-point wavefunction obtained by the PWFRG. By sweeping the 
field rj, we obtain a p — rj curve which plays the role of the "equation of state" . Actual 
calculation is very similar to that of the magnetization curve for spin chains f22| , |24}| (see also 

!)■ 



The Ising magnetization M which is calculated by < 07, > where < • > denotes the 
thermal average. The adsorption coverage B is connected to M by = (M + l)/2. 

III. ADSORPTION MEDIATED INTER-STEP ATTRACTION 



In the RSOS-Ising coupled system, we find, at low temperature, the crystal shape shows 
the "first-order" transition at the facet edge, that is, the surface gradient p jumps from 
to a finite value p c (T, H) at the facet edge. In Fig. 2, we show p — rj curve and M — rj 
curve calculated by PWFRG method at J/e — 0.15, H/e = —0.01. Both M and p jump at 
facet edge at low temperatures. 

We also have made Monte Carlo calculation at J/e = 0.15, H/e = —0.01 and temperature 
ksT/e = 0.4. In Fig. 3, we show the top view of RSOS surface (above) and the adsorbates 
configuration (below). Steps get together to form a "quasi-facet" , and the surface gradient 
of the quasi-facet agrees with the p c (0.4, —0.01) in Fig. 2. The vicinal surface of a mean 
surface gradient p breaks up into two regions: one with p = and the other with p = p c . This 
phenomenon is precisely the phase separation characteristics to the first-order transition P? . 
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We should remark that the quasi-facet with p = p c is not a "true" facet in the ECS drawn 
by the Wulff construction. The facet must be the singular surface corresponding to the cusp 
orientation of the Wulff figure p8| , |26|| ; this is not the case for the quasi-facet. Although the 
quasi-facet is rough, the boundary between the p = surface and the quasi-facet ("quasi- 
facet edge") is "sharp". 

From a microscopic point of view, the above thermal step bunching implies the existence 
of attractive interaction between the steps. To confirm this, we calculate a "terrace width 
distribution function" P{1) of the RSOS-Ising system (number of steps = 2) by Monte Carlo 
method up to 4 x 10 7 MCS (Fig. 4). We denote locations of step ledge at y of the step 1 
and the step 2 by Xi(y) and x 2 (y), respectively. The linear size of the square lattice in the x 
direction and the y direction is denoted by L x and L y , respectively. We impose the periodic 
boundary conditions along both x- and y-directions. To be precise, the "periodic boundary 
condition" along the x-direction is the one where the uppermost terrace is connected to the 
lowermost terrace. The "terrace width" I is calculated to be |x 2 (y)— Xi(y)\ or \xi(y)+L x — x%\. 

The solid curve in Fig. 4 clearly shows double peaks with very narrow peak width (this 
shape is unchanged under the increase of L y ). For J = (broken line) where only the 
entropic repulsion exists between the steps, the P{1) drops around I ~ or / ~ L x , with no 
peaks. It is clear, therefore, that the attraction exists between steps for J ^ case. 

Origin of this attractive force is the fluctuation of the Ising spins (number density fluctu- 
ation of the adsorbates). After coarse graining, Ising system can be described in terms of the 
4 -like continuous "massive" field [2{| except for at criticality (the transition temperature of 



the Ising system at zero field is 0.15 x 2/ ln(l + \^2) ~ 0.3404). Due to the coupling term of 
the form —aa\Ah\, there appears effective force between the height differences {Ah}. The 
force mediated by the massive field is short-ranged, with the range being ~ £ (£: correlation 
length of the Ising model). In other words, the adsorbates play the role of photons (for 
Coulomb force) or Yukawa mesons (for nuclear force). This attractive interaction leads to 
the instability against step bunching. 

At the end of the first-order transition line, there usually appears a critical point where 



the second order transition occurs. In Ref. |L9] We find that B becomes zero at a temperature 
where the surface free energy has the form 

f(p) = f(o) + 7 p + CfP + (higher order). (3.1) 

The GMPT behavior ( |1.1| ) breaks down at this point. In jB|, we have presented two possible 
mechanisms for the vanishing of B ("stiffening transition"): (1) "single-particle" mechanism, 
and (2) intrinsically many-body mechanism. From the present study, we conclude that the 
second mechanism is the case, and that the "stiffening transition" is naturally interpreted 
as the one which is very much like the critical point in the liquid-gas transition. 
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FIGURES 

FIG. 1. (a) The vertex of the RSOS model. A vertex is drawn along the boundary between 
four adjacent SOS columns. At the legs of the vertex, variables [i, %, 5 and v are assigned as 
hi — hi, h± — hs, hi — h% and hi — h±, respectively. These variables take {1,0,-1} and should 
satisfy n + x = S + u. (b) The vertex representation of the Ising-RSOS coupled system. Ising spins 
which take ±1 are described by open circles on the legs of vertex. 

FIG. 2. The PWFRG results for the p - (3r) curves with a = 0.5, H/e = -0.01 and J = 0.15 
where (3 = 1/k^T. The temperature of each curve is k^T je = 0.3,0.4, and 0.5 from right to left. 

FIG. 3. The snapshot of the Monte Carlo simulation at 10 6 MCS. Above: Configuration of 
RSOS vicinal surface; a black dot shows Ah = —1, a white dot shows Ah = 1, and gray region 
corresponds to Ah = 0. Below: Configuration of Ising system (adsorption layer); a white dot shows 
(7 = 1 (occupied) and a black dot shows a = — 1 (empty), a = 0.5. H/e = —0.01. k^T je = 0.4. 
Size: 408 x 120. We have 24 steps on the surface. 

FIG. 4. The "terrace width distribution" calculated by the Monte Carlo method averaged 
over 4 x 10 7 MCS. a = 0.5. H/e = -0.01. k B T/e = 0.4. Size: 70 x 120. The bold line and the 
broken line shows J = 0.15 and J = 0, respectively. 
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Figure 1 N. Akutsu, Y. Akutsu and T. Yamamoto 
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